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Abstract 

The analytic properties of scattering amplitudes provide important information. Besides the cuts, the poles and zeros on the 
different Riemann sheets determine the global behavior of the amplitude on the physical axis. Pole positions and residues allow 
'for a parameterization of resonances in a well-defined way, free of assumptions for the background and energy dependence of 
the resonance part. This is a necessary condition to relate resonance contributions in different reactions. In the present study, we 
determine the pole structure of pion-nucleon scattering in an analytic model based on meson exchange. For this, the sheet structure 
of the amplitude is determined. To show the precision of the resonance extraction and discuss phenomena such as resonance 
interference, we discuss the S \ \ amplitude in greater detail. 

Key words: meson-baryon scattering, baryon resonances, analytic continuation 
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X. Introduction 

The global properties of a scattering amplitude are deter- 
mined by the kinematics of the reaction, leading to branch cuts 
associated with the opening of reaction channels of stable or 
unstable particles. The thresholds of two particles or quasi-two 
particles being on-shell are characterized by the branch points. 
There is a righthand and a lefthand cut associated with s chan- 
nel and crossed channel processes, and there can be resonances, 
bound and virtual states. Resonances and virtual states are asso- 
ciated with poles on unphysical sheets. Thus, an analytic con- 
tinuation along the various branch cuts is mandatory to access 
the resonance poles. 

, The nN — > jiN transition is one of the most precisely mea- 
sured reactions. It provides detailed information about the 
baryon spectrum, which is presently under experimental inves- 
tigation, see e.g. Ref. 111]. Masses, widths, and decays of bary- 
onic resonances allow for tests of models of the internal struc- 
tures of the nucleon and its excited states. 

Yet, the resonances in the second resonance region have to 
be disentangled. Thus, most of the prominent resonances listed 
by the PDG have been obtained by partial wave analy- 
ses Jim in followed by a model dependent analysis of the par- 
tial wave amplitudes e.g. in terms of a background and Breit- 
Wigner resonances JHIa]. In the energy range between 2 and 
3 GeV, presently under experimental investigation, resonances 
start to overlap and the background may show some non-trivial 



structures. This situation calls for more sophisticated theoreti- 
cal analyses as e.g. carried out in the partial wave analyses of 
Refs. ySllTI] where poles in the complex plane of the scattering 
energy are determined. 

In the case of overlapping resonances and resonances near 
thresholds, partial decay widths can only be extracted within 
models. When different reactions are analyzed simultaneously, 
this procedure becomes questionable. We therefore stress that 
the only sensible parameters that encode the resonance proper- 
ties are the various pole positions and residues. 
Models of the K matri x type | 



meson exchange models 




IJand 



unitary amplitudes that have been constructed in the past to ac- 
cess pion-nucleon scattering. In unitarized chiral perturbation 
theory, resonances can be generated by the non-perturbative in- 
teraction of mesons and baryons without the need to e xplicitly 
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In this context, one has to identify observables that allow to 
distinguish between hadronic molecules and elementary states. 
Ref. 13011 shows how to incorporate explicit resonance fields 
in a unitarized extension of chiral meson-baryon dynamics and 
stresses the importance of an accurate quantitative reproduction 
of the experimental data which is prerequisite for a determi- 
nation of the pole position in the complex plane. Weinberg's 
method to decide on the molecular nature of bound states 1371] 
which is valid for i-wave states close to threshold has been gen- 
eralized to resonances in Ref. llssll . 

In some of the approaches to meson-baryon scattering, the 
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amplitude has been analytically continued to the complex plane 
to extract the pole positions and parameters, see e.g Refs. 

|32 [ [33 ! [35 ll and the recent work of Ref. ll39[l within the me- 
son exchange framework. In the present study, we extend the 
amplitude of the Jiilich model to the various Riemann sheets in 
the complex plane of the scattering energy i'^^ = z- 

The Jiilich model is an analytic coupled channel model based 
on meson exchange that respects two-body unitarity. This 
model has been developed over the past few years 11911 . with 
its current form, as used in this study, given in Ref. [2Q|. For 
the convenience of the reader, we point out the main ideas in 
the following. 

The coupled channel scattering equation IS 
solved in the JLS basis, given by 

{L'S'k'\T'J^.\LSk) = {L'S'k'\V'J,\LSk) 
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where J {L) is the total angular (orbital angular) momen- 
tum, S (/) is the total spin (isospin), k{k' , k") are the incom- 
ing(outgoing, intermediate) momenta, and /y, v, y are channel 
indices. The incoming and outgoing momenta can be on- or 
off-shell. The integral term involves a sum over all interme- 
diate possible quantum numbers and channels contained in the 
model. 

The integral term can be abbreviated as VGT where G is the 
intermediate meson-baryon propagator of the channels with sta- 
ble particles jiN and rjN, given by the fraction in Eq. ([T]i. For the 
channels involving quasiparticles, crN, pN, and ttA, the propa- 
gator is slightly more complicated il8lll9tl . The pseudopoten- 



tial V iterated in Eq. ([T]) is constructed from an effective inter- 
action based on the Lagrangians of Wess and Zumino 14 ll 14211 . 
supplemented by additional terms lfl9[ Eoll for including the A 
isobar, the w, meson, and the cr. The exchange potentials 
V are partial wave projected to the J, L, S basis, where Eq. ([T]) 
is solved. Note the potentials V and amplitude T appear half- 
oflF-shell in the integral term of Eq. ([T]i. 

The exchange potentials mentioned above contribute to the 
non-pole part. The pole part is given by baryonic resonances 
up to y = 3/2 that have been included in V as bare s channel 
propagators. The resonances obtain their width from the rescat- 
tering provided by Eq. ([T]i. 

In this study, the partial wave amplitudes T'Jj^ will be an- 
alytically continued, rather than the full amplitude or invariant 
amplitudes A-(s, u), B^{s, u). Partial wave amplitude have a 
rather involved structure below the nN threshold which is dis- 
cussed in Appendix lAl While the continuation for the channels 
with stable particles is straightforward, the effective nnN chan- 
nels require special attention. It is known ^431 ] that the quasi- 
two particle singularities induce novel structures in the ampli- 
tude, i.e. additional branch points in the complex plane, apart 
from the righthand nnN cut along the physical axis. The result- 
ing sheet structure is non-trivial and should be fully taken into 



account; these additional branch points induce large variations 
of the amplitude in their surroundings and have a large impact 
on pole positions and residues. For the Roper channel, the role 
of additional sheets has been discussed in Ref. 144 1. 

In some studies, the properties of channels with unstable par- 
ticles have been modeled by folding the invariant mass distribu- 
tion of the unstable particle with the stable particle propagator 
as e.g. in Refs. 11451 14611 . We have developed related techniques 
based on a Lehmann representation for the righthand cut. While 
the discontinuity can be expressed in this form, an explicit proof 
of analyticity of the obtained continuation is difficult within our 
framework. 

For the Jiilich model, the various sheets are accessed through 
contour deformation of the momentum integration. A similar 
method has been recently applied for another model of the me- 
son exchange type . 

Another issue of relevance in this context is the question, if 
it is possible to remove the hadronic contributions from observ- 
ables in a model independent way to allow access to quanti- 
ties that can be identified with those calculated from the quark 
model 1 47 , 48 1 — see also Ref. ll49ll for a recent discussion of 
the subject. Such an analysis assumes that a clean cut separa- 
tion of pole and non-pole parts is possible. 

In a recent study ifsoll this question has been addressed. 
The conclusion, drawn from the behavior of the P33 partial 
wave with the A(1232) resonance, was that the residues pro- 
vide a meaningful expansion parameter of the resonance am- 
plitude, while dressed and bare vertices depend on the regu- 
larization scheme and the model dependent decomposition of 
the amplitude into pole and non-pole contribution according to 
T - + T^^. This separation is widely used in the literature, 
see e.g. Refs. ifsil 15211 . Second, the non-pole part itself can 
contain poles, that interact with the pole part in a non-trivial 
way, as discussed in Ref. ifsoll . In the present study, we test the 
Di3 and other partial waves in the light of these questions and 
come to similar conclusions. 

After all, especially for the comparison of different exper- 
iments, the poles and residues of the S-matrix are the relevant 
quantities. The Jiilich model is derived within a field theoretical 
approach from Lagrangians obeying chiral constraints. It pro- 
vides a precise description of the partial wave amplitudes. Fur- 
thermore, other than in most K matrix approaches, analyticity 
is respected. These ingredients qualify for a reliable extraction 
of resonance properties. 

This paper is organized as follows. In Sec. 12.11 the well- 
known analytic structure of the propagator of two stable parti- 
cles is discussed. The more complicated properties of quasi- 
twoparticle intermediate states are discussed in Sec. 12.21 Re- 
sults for pole positions, residues, zeros of the amplitude and 
transition strengths are provided in Sec. [3] 

2. Analytic structure of the scattering amplitude 

2.1. The propagator of two stable particles 

The analytic continuation of the scattering amplitude will be 
discussed in several steps. In this section, the well-known struc- 
ture for the propagator of two stable particles is presented. The 
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Figure 1 : Real (dashed) and imaginary (dotted) parts of the cr self-energy at 
= MeV. 



connection between analytic continuation and contour defor- 
mation of the momentum integration will be pointed out. The 
contour deformation is the key to the analytic continuation of 
the effective nnN channels crN, pN, and nA, discussed in Sec. 



To simplify the discussion, we choose pion-pion scattering 
in the L - I - Q channel. Here, the scattering problem is de- 
scribed as a one-channel amplitude with a propagator of two 
stable pions and one explicit resonance, the cr(600) (7^ = 0"^) 
meson. In the Jiilich model, this amplitude appears in the con- 
struction of the crN propagator and also in the cr f-channel ex- 
change. The example will serve to demonstrate the analytic 
properties of the propagator of two stable particles. The follow- 
ing discussion applies qualitatively also to the ttN and 7]N chan- 
nels of the Jiilich model, as well as to other tttt and ttN channels, 
namely the p(770) in the pN propagator and the A(1232) in the 
ttA propagator 

The cr can be described with a dressed TOPT (time-ordered 
perturbation theory) propagator according to 
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(2) 



where = 900 MeV is the bare cr mass, k = \k\ is the cm. 
momentum, i.e., A: = in the cr cm. frame, z' is the energy 
boosted to the nn cm. frame (z' = z at k - is the only case 
needed in this section), and 11^ is the cr self-energy given by the 
loop 
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where v'^'"'{q,k) is given in Eq. (3.10) of Ref. 11 1811 . For the p 
and A self energies, the vertices are given by Eq. (20) of Ref. 
Ill and in Ref. respectively. 

The cr self-energy for real energies z is shown in Fig. [T] 
The imaginary part becomes finite and negative at the two-pion 
threshold while the real part shows a cusp at this point. To un- 
derstand the analytic properties of the cr self-energy, consider 
complex values of z in Eq. (O, evaluating the integral along 
a straight path from zero to infinity along the real q axis. The 
self-energy as a complex function of z is shown in the first Une 
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Figure 2: The two Riemann sheets of the cr self-energy [arb. units] as a function 
oiz [MeV]. The left column shows the real and the right column the imaginary 
part of Flo-- 



of Fig. |2] The result is called the first sheet of the self-energy. 
There is a righthand cut extending along the positive real z axis 
starting at the energy where both pions can go on-shell, i.e. 
■Zthresh - 2 m„. The cut is at z - ie, i.e. Ho- on the real z axis 
is analytically connected to Ho- in the upper half plane. 

Along the righthand cut, the amplitude can be analytically 
continued. For this, the standard procedure is to calculate the 
imaginary part at the real axis by evaluating Eq. ((SJ with the 6 
function 6{z - 2Eq) resulting in 



T T-r ■"'/on^on^on 2/ > ;\ 
ImHo- = V (^on' ^) 



(4) 



where i^on is the on-shell three-momentum in the cm. frame of 
the stable particles 1 and 2 (in this case two pions) and E^a is 
the on-shell energy of the particles. As the quantity 



?on = TT V(z- - (nil - '«2)^)(z^ - im + mf) 
2z 



(5) 



is two-valued itself, we need to distinguish the two Riemann 
sheets of ^on uniquely according to 



^on 



if Im qo 
else 



< 



(6) 



With this definition, q^^ is real and positive on the real axis 
above threshold; it has a cut along the real axis at z - ie above 
Zthresh, whcrc the branch point is located, and is analytic other- 
wise. 

The analytic continuation to the second sheet is given by 
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Figure 3: First (dashed lines) and second sheet (solid lines) of the cr self-energy 
Ilo- [arb. units]. The upper row shows a slide for Rej above the nn threshold 
(note the discontinuity at Im z = 0), the middle row for Rej = 2m,r and the 
lower row for Re z below the threshold. 



i.e., by subtracting twice the imaginary part of Eq. (HJi. In Fig. [3] 
we show the two sheets hJ^^' s n,j- (soUd Unes) and IlJ^' (dashed 
lines) as a function of Imz for fixed Re (z). In Fig. |2]the same 
situation is shown in a three-dimensional plot. 

The two sheets are analytically connected along the right- 
hand cut. Note that below threshold the second sheet is not di- 
rectly connected to the physical axis but only via paths around 
the branch point at Zthresh- Above threshold, there is a direct 
connection from the second sheet in the lower half plane to the 
physical axis. 

In the following we sketch the proof of analyticity of the pre- 
scription from Eq. (|7]i. Although this proof is trivial it will 
help us understand the procedure of the analytic continuation 
for unstable particles in Sec. 12.21 

Consider the pole of the integrand from Eq. (O with respect 
to q in the right q half plane. Its position, for different z values, 
is indicated in Fig. |4]with crosses. For z values with positive 
imaginary part -shown as case (a)- the pole is in the upper q 
plane and the q integration can be carried out on a straight path 
from zero to infinity. For real z -this is the relevant case to 
evaluate observables- the pole is on the real axis [cf. case (b)] 
and one can deform the integration contour with a half circle as 
shown in Fig. H] The evaluation of the infinitesimal half circle 
returns the imaginary part of the self-energy. Once Imz < 0, 
the pole is in the lower half plane. This crossing of the pole 
over the integration contour corresponds to the appearance of 
the right-hand cut. Yet, one can obtain an analytic continuation 
for negative Imz by defining a deformed path shown as case (c) 
in Fig. m 

From the construction of the path in case (c) it becomes clear 
that this continuation is indeed analytic. From the form of the 
path it is also clear that the result of the integration corresponds 
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Figure 4: Integration contours (dashed lines) for the loop of two stable particles 
imr). For different the poles are marked with a cross. The arrows indicate the 
changes of the pole positions for the cases Im z > 0, = 0, < 0. 



to an integration along a straight path on the real axis (corre- 
sponding to the first sheet), plus the integral taken on the full 
infinitesimal circle, that returns the residue of the integration. 
The contributions from the vertical parts of the path (c) in Fig. 
|4]cancel. The final result is then exactly given by Eq. as an 
explicit evaluation of the residue shows. 

To finish this discussion, let us again consider the poles in 
the ^-plane. The pole at qp in the right half plane is simple; it 
is accompanied by another simple pole at -qp. However, for 
the energy Zh, - m\ + mi there is only one double pole in the 
momentum plane, situated at ^ = 0. The energy Zh^ , however, 
is threshold and branch point of the Riemann surface. 

Concluding, a branch point for integrals of the type of Eq. 
(O is induced whenever the two poles in the ^-plane coincide 
and form a double pole. We will also encounter this situation 
in the more complicated case of the effective nnN propaga- 
tors, where the classification of the Riemann surface in terms 
of branch points is more complicated. 

Eq. ^ gives the analytic continuation of the self-energy of 
two stable pions. In the Jiilich model, there are the two channels 
with stable particles jiN and rfN. The analytic continuation of 
the amplitude with respect to a channel (mn) given by a stable 
meson m and a stable baryon n is closely related to the second 
term in Eq. (|7]i, i.e. the discontinuity along the righthand cut. In 
particular, we obtain the scattering equation for the amplitude 
r*^' on the second sheet. 



{qcd\T^^' - VVlab) 
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6G = 



2mq>^(mn)El-Ef, 



X {qcd\V\qUmn)}{q^^^(mn)\T''%,i,} (8) 



with £,„„ - E„i + E„ and q^^ from Eq. (|6]l. Indices of quantum 
numbers and the angle integration have been suppressed in Eq. 
(HJl. Eq. ^ has a similar form as the scattering equation on 
the first sheet from Eq. ([TJ, except for the additional term SG, 
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which is on the real axis given by the discontinuity of the two 
particle propagator l/(z - E,„„). In Eq. ([8]), the indices (ab), 
(cd), and (mn) indicate the incoming, outgoing, and intermedi- 
ate channels, respectively. 

Note that the matrix elements V and T*^' in the expression 
for 6G appear in on-shell kinematics which allows for an imple- 
mentation of the term at the on-shell point of the Haftel-Tabakin 
scheme 15311 that is used to solve the scattering equation. 

One can now also search for poles of the cr(600), p(770), and 
A(1232) as they appear in the effective ttttN propagators of the 
Jiilich model. Simple but analytic models have been developed 
to obtain a good empirical parameterization of the tttt and nN 
amplitude which are needed as input just for a parameterization 
of the effective nnN channels. This purpose is well fulfilled by 
these models. 

For the pole search, we analytically continue the propagator 
from Eq. (|2]i, with n[^* from Eq. according to 
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(9) 



the (T, p, and A to allow for a realistic description of the effec- 
tive nnN channels as they appear in nN scattering, the principal 
objective of this study. 

2.2. Propagator with unstable particles 

The analytic continuation for the effective ttttN channels crN, 
pN, and ttA is different from the channels nN and rjN discussed 
in the previous section. There, we have seen that a deformation 
of the integration contour as shown in Fig. |4] leads to an ana- 
lytic continuation. In practical terms, one can simply add the 
discontinuity using the 6 function which is equivalent. For the 
effective nnN channels, the key to the analytic continuation is 
the contour deformation. 

For the discussion, we focus on the crN channel. The pN and 
nA channels can be treated analogously as discussed at the end 
of Sec. 12.41 The additional complication that arises here is that 
the unstable particles appear in a moving frame. Within TOPT, 
the three dimensional formalism used here, the cr self-energy 
evaluated for a finite sigma three momentum k reads 



as A: = in the cr cm. frame. A pole on the second sheet of C C dx {v°''"'{q k)Y 
the ;r;r scattering amplitude corresponds to a pole of G[^'(Zo-). t\a-{z,k) - I q'dq I — — — —, (10) 
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We find the pole of the cr at z^J. = 875 - 232 / MeV. The cr pole 
is situated relatively high in energy compared to the high pre- 
cision determination in Ref. Q of z° = 441 - 272 / MeV. At 
this point the model could be improved by developing a scheme 
with derivative coupling that allows for a cr pole at lower ener- 
gies. In any case, a satisfying parameterization of the nn scat- 
tering amplitude in the cr channel is obtained. In Fig. 10 of Ref. 
lfl9ll the quality of the fit is shown. There are some minor devi- 
ations, probably due to the displaced cr pole, but for the present 
purpose, this accuracy is sufficient to provide a realistic descrip- 
tion of the three body phase space, starting at z = 2m„ + up 
to the maximum energies considered of around z ~ 1.9 GeV 
In the complex plane, a change in the pole position of the cr 
would lead to a change of the position of the branch point (cf. 
Eq. (flSTl): however, as the pole is quite far in the complex 
plane, no major changes from a modified cr pole are expected 
in the region where poles of baryonic resonances are searched 
for(r/2 < 150 /MeV). 

For the p, a similar analysis can be made. There are no poles 
on the first sheet and one pole on the second sheet at Zp = 763 - 
64 ; MeV. This corresponds to ap width of F = 128 MeV which 
is slightly smaller than the standard value of F = 150 MeV 

For the A, the pole Hes at z° = 121 1 - 37 ; MeV, i.e. the A has 
a width of F = 74 MeV which is smaller than standard values of 
1 10-120 MeV H. The pole position of the A(1232) in the nA 
propagator is not identical to the position of the A(1232) in the 
nN i-channel exchange, that will be determined later (cf. Ref. 
IS or Table EJ. The consistency of the model could be im- 
proved at this point, i.e. consistent pole position for both cases. 
However, the discrepancies are minor. As we have stressed be- 
fore, at this point a high precision fit of the amplitude is not 
required; a reliable parameterization of the input for the effec- 
tive nnN channels is sufficient. 

Thus, the nn and nN phase shifts are sufficiently well de- 
scribed (cf. Fig. 7 of Ref. ifli]. Figs. 9 and 10 of Ref. Ill) by 



with tj^ - ■\j{q + k/2)~ + ml and x for the angle enclosed by q 

and k. This expression reduces to the one given in Eq. (O for 
k = 0. Note, fla- is a function of only the modulus of k and not 
its direction as a result of the x-integration. It is instructive to 
expand the denominator of the self-energy for small values of 
k. We then get 



■2u,-k-l{Au,)+0{k^lajl), 



with cl>„ — -\Jq~ + ml. Thus, through the boost momentum k the 
energy available for the cr gets reduced by A^/(4<y), the kinetic 
energy of the two pion system. We may therefore parameterize 
the crA^ propagator as 118 1 



gcrNiz, k) 



z- ^ml + Ar2 - ^J{m'^,)^ + k^-UAzAz,k),k) 



Go-Niz) = J dkk"- F(k)ga-N(.z,k), 





zAz,k) ^z + ml - ^k^ + (m0)2 - ^/F7i 



(11) 



with the self energy Ho- from Eq. (O; m° is the bare cr mass and 
F is a regulator we introduce for the discussion of this section 
and which is absent in the Jiilich model [cf. Sec. 12.511 . The term 
guN is given by the TOPT crA^ propagator that includes the self 
energy of the cr meson. The energy parameter Zo- implies that 
the energy available for the nn subsystem is reduced compared 
to the total energy not only by the kinetic energy of the cr, but 
also by the energy of the nucleon propagating simultaneously. 



5 



9r^jKl 



(7 



I 
I 



N 



Figure 5: Example of an intermediate (tN propagator in the rescattering pro- 
vided by tlie scatteiing equation 0. 
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Figure 6: Real (dashed) and imaginary (dotted) parts of the crN propagator Ga-N 
for real energies z [arb. units]. 
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Figure 7: In the k plane (upper figure), the path of the ^-integration (long dashed 
line) can be deformed (short dashed line). In the Za- plane, this con'esponds to 
dift'erent paths for different values of z (lower figure). In the Ztr plane, there is 
a branch point at Zo- = 2m^. The corresponding cut can be chosen in different 
directions, e.g. (a) or (b). 



We checked numerically that treating Ho- as in Eq. (O vs. keep- 
ing the full dependence of 11^ of Eq. ( fTOl i has only a very small 
impact on observables. Since using Eq. (fTTT i saves an angular 
integration we will use this expression in what follows. 

In Eq. ( fTTT i. we have explicitly included the integration over 
k, the loop momentum of the crN loop. As an example of such a 
loop, we show in Fig. |5]two pion exchanges. The exchange pro- 
cesses are k and angle dependent; they induce three-body cuts 
and additional analytic structures, on top of the structure given 
by gijj^. This issue will be further discussed in Sec. I2.5l where 
also the implementation of the analytically continued propaga- 
tor into the full Jiilich model is given. See also AppendixlAlfor 
a discussion of the short nucleon, circular and other cuts. 

In the following, the analytic structure of G,jn as defined in 
Eq. (fTTT i is determined. For real energies z, the propagator G,jn 
is shown in Fig. |6] The real and imaginary parts can be com- 
pared to the case of stable particles as shown in Fig. [T] One 
observes similar structures, which are, however, smeared in en- 
ergy due to the finite width of the cr, or equivalently, a finite 
in Eq. (fTTT i. In particular, the influence of the "threshold" at 
z ~ m% + niN ~ 1.9 GeV is still visible and reminds one of the 
corresponding structure in case of stable particles. 

The 7T7T self-energy 11^^ in Eq. ( fTTT i has a well-known right- 
hand cut along the real Za- axis as derived in Sec. 12.11 It is clear 
that this induces also a cut in the full propagator Go-n- 

In the case of the propagator of two pions the discontinuity 
arises from the position of the pole with respect to the integra- 
tion contour, as discussed following Fig. H] In particular, for 
real z, the pole in the ^-integration lies on the real axis. In the 
present case, the pole lies, for real z, far in the complex plane 
[c.f. Eq. (fTTTil. and the cut of Go-n on the real z axis, starting at 
z = 2m„ + niN, is entirely induced by the cut of IIo- itself. 



Thus, for an analytic continuation of Go-n, one first has to 
analytically continue Ho-. For this, we consider the k integration 
in Eq. (fTTT i. The integration contour in the complex k plane 
can be deformed as long as the limits (0, oo) are unchanged. 
Second, the deformation must be along analytic regions of the 
integrand and must not cross poles, branch points, or cuts. This 
guarantees that the result of the integration along the deformed 
contour is unchanged. 

To ensure these conditions, we have to check which are the 
Zo- regions in Do- that correspond to the contour deformation in 
the k plane. A path in the k plane corresponds to a path in the 
Zcr plane as given by the transformation Zcr(z, k) from Eq. (fTTT i. 
We show in Fig. [T] schematically several deformed paths in the 
k- and Zo- -plane. 

First, consider the undeformed path for a value of Imz = 
Im ZiT > 0, i.e. a straight path in the k plane from zero to infinity. 
This is indicated with the long dashed lines in Fig. [T] The 
integral along the deformed path (short dashed lines), given that 
the endpoint is the same as for the undeformed path, will return 
the same result because the deformation is over a region in Zo- 
where Ho- is analytic. 

Next, consider the case Imz = Im Zo- < 0. The righthand cut 
of Ho- is along the real Zo- axes, indicated as direction (a) in Fig. 
[T] The righthand cut begins at the branch point Zo- = 2m„. Thus, 
the integration along the undeformed path induces a cut in the 
full Ga-N once Imz changes sign. 

The analytic continuation of G^n along this cut, for Imz < 0, 
is obtained as follows: If Imz changes from positive to negative 
values, the integration contour is deformed as indicated with 
the short dashed line in the Zo- plane in Fig. |7](case Imz < 0). 
However, this is possible only if simultaneously the cut of 11,^ is 
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Figure 8: Imaginary part of the self-energy n*f' [arb. units] as a function of 
Za- [MeV]. Tlie cut is put in the negative Imio- direction. With this choice, the 
self-energy is analytic on the real axis above the branch point. 



moved from direction (a) to direction (b). This ensures that the 
integrand is always analytic. The continuation obtained through 
this contour deformation is the unique analytic continuation of 

GcrN- 

The change of direction of the cut from (a) to (b) redefines 
the self-energy according to 



(2) 



if Im Zo 
else 



< and Re Zo- > 2m„ 



(12) 



where Ho- and 11^' are given in Eqs. (O and (|7]l, respectively. 
The self-energy 11^^ is shown in Fig. [8] Indeed, the cut is now 
along direction (b) and 11**' is analytic on the real axis above 
the branch point. 

Second, as shown above, the integration path has to be de- 
formed; for simplicity, we have decomposed the path into 
straight pieces in the k plane in such a way that one interme- 
diate edge point corresponds to Zcrikmd = 2 m^r + 100 ; MeV. 
The start and end points are always given by ^ = 0, ^ = oo. 

For Re z < 2m„ + we obtain Re Za < '^mn and the starting 
point of the integration at A: = lies below the cut at position 
(b). Yet, one simply rotates the cut direction further than -90° 
and can apply the method as before. 

The analytic continuation of Go-n to the lower z half plane, 

(21 

called G^f^, has been evaluated. A quantitative definition of 
G^^ will be given below [cf. Eq. (fT6Tl1. To continue G^^ to the 
upper plane, one can adapt the method of contour deformation, 
or simply utilize the known continuation in the lower half plane 
and take advantage of the general analytic property 



(13) 



that holds also for the second sheet. 

In Fig. |9]we show the analytic continuation as slices along 
the Im z direction at fixed Re z- The solid lines show the first 
sheet G^j*^ for Imz > and G[^^ for Imz < 0. As Fig. |9] shows, 
the transition between both sheets is indeed analytic at Im z = 0. 
The dashed lines show, for Imz < 0, the first sheet G[^'^. The cut 
in G[^^, situated at Imz - 0, appears in Fig. |9]as discontinuity 
for the imaginary part and a non-analyticity ("cusp") for the real 
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Figure 9: Real (left) and imaginary (right) part of Go-w [arb. units]. The plots 
show slices along the imaginary z direction for Re z fixed at 1.4, 1.75, 1.87, and 
2.1 GeV (top to bottom). The solid lines show the analytic continuation to the 

a-N 
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part. The second sheet G)~^ for Imz > 0, obtained through Eq. 
(fTSl) . is not plotted in Fig. |5] 

2.3. Additional branch points in the complex plane 

The analytic continuation of Im Go-n is shown with the solid 
lines in the right column of Fig. |9] While, at Re z = 1 .75 GeV, 
the imaginary part shows a sharp rise as Imz < -250 MeV, 
only 120 MeV above at Rez = 1.87 GeV the continuation falls 
rapidly beyond Imz < -250 MeV. The real parts show also 
rapid changes in this z region. This is a sign that there is an 
additional structure. In the following we show that the new 
structure is induced by a branch point. Such additional branch 
points are known since long 11431 14411 . 

In Sec. 12. II it has been shown for the case of propagators 
of stable particles that branch points are induced whenever two 
poles in the complex plane of the momentum integration over 
q coincide, i.e., the denominator has a double zero in the mo- 
mentum plane and ^ = 0. For the stable case, this situation was 
given at threshold. In the present case, we have to inspect the 
denominator of Eq. (fTTl i and search for poles with respect to 
k. This is, in general, only possible numerically. However, we 
know that for a branch point, k -0. Furthermore, the condition 
^ = and the fact that we consider the analytic continuation 
of G[^^ imply that the cr self-energy is evaluated on the sec- 
ond sheet, n[^^ according to Eq. (fTST i [c.f. Fig. [Tj. Then, the 
denominator of Eq. ( fTTT i has the double zero at 



Im z 
▲ 



Z - niN - nil - n<?'(2 - mN, 0) = 0. 



(14) 



On the other hand, the pole of the cr itself is given by Eq. (|9]i; 
using that Zo- - z — niff at ^ = 0, we obtain for the position of 
the branch point, called b2 in the following. 



Zb2 



Zn + niN 



(15) 



where Zp = 875 - 232 / MeV is the pole position of the cr in the 
complex Za- plane. 

Thus, branch points in the complex z plane of the effective 
nnN propagators are directly related to the pole of the unstable 
particle. An unstable particle cr induces branch points in the 
crN propagator. These branch points are always on the second 
sheet of Go-n, as we have seen in the derivation of Eq. ( fTSl l. 
because this is where the resonance poles are. Furthermore, as 
there is only one pole on the second sheet of the cr propagator, 
there are no further induced branch points of G,rN in the z plane, 
apart from b2 [and fe^ in the upper z half plane due to Eq. (fT3]l1. 

For the pN and ttA propagators, relation (fTsl l holds as well, 
with the corresponding masses and pole positions; in all cases, 
the validity has been confirmed numerically. 

Relation ( fTSI l has been derived for the specific form of the 
propagator from Eq. (fTTT i. Yet, the existence of the branch 
points bo, b'2 follows in general from the existence of the pole 
in the scattering amplitude of the unstable particle. Second, the 
validity of Eq. (fTsl l does not depend on the special form of the 
transformation Zo- from Eq. (fTTl i: the condition Za- - z - 
at A: = that led to Eq. ( fTTI ) simply reflects the opening of the 
jtttN threshold at z = 2rn„ + mt^ and is obeyed in general Q. 
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Figure 10: The full analytic structure of the crN propagator Go-n- The branch 
point b 1 is at the nnN threshold and connects first and second sheet. The branch 
points bi and i)^ are located at zbi = Zp + '«« (Zh' = zj, + and connect 
second with third, and second with fourth sheet, respectively. The lines i | and 
S2 indicate slices plotted in Fig. [12] .S3 indicates the viewpoint of the plot in 

Fig.\n\ 



The resulting analytic structure is shown in Fig. \TQ\ There is 
a branch point bi, located at zt^ = 2ot^ + m^, which connects 
first and second sheet 0, both of them with a cut along the real 
axis. As we have seen previously, this branch point and its cut 
are induced by the cut of the cr self-energy Ilo-. The additional 
branch points b2 and lie in the complex plane, both of them 
on the second sheet, and induce the two additional sheets three 
and four. 

In Fig. [TT]we show the real part of Go-n around b2 from a 
viewpoint 53 as indicated in Fig. [TO] There is an intersection 
of real parts visible, starting at ^2 and extending towards pos- 
itive Re z values. For a propagator of stable particles, the real 
parts of the two sheets intersect along the real axis, or righthand 
cut, as can be seen in the upper left plot of Fig. [3] It is, thus, 
straightforward to define the cut belonging to b2 in the same 
way, i.e. at the intersection of the real parts. This is indicated 
in Fig. [TO]with the red curved line. 

In Fig. [12] the four Riemann sheets along the two slices si 
and S2, as indicated in Fig. [TOl are shown. The slices are lo- 
cated slightly below and above b2. In order to understand the 
structure of the branch point, one can follow paths around it. 
For example, coming from the real z axis on the second sheet, 
one can follow si and pass by b2 below (thick dashed line in 
Fig. [T2I 1. Then, one can move to S2 and move back, towards 
the real axis, passing b2 above (thin dashed line). As Fig. [12] 
shows, one is then not any more on the second sheet, but on 
sheet 3. Alternatively, one can start at the real z axis on sheet 2, 
along S2 (thick black solid line), move around b2 and return on 
si (thin solid line), and also get to sheet 3. Note the presence of 
a complex conjugate structure that follows from Eq. ( [T3] ), as- 



'The existence of bj and its general properties are not discussed in Ref. 

^The counting of the sheets refers in this section to the one channel case of 
o-N. 
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Figure 11: Real part of Ga-N [ai'b- units] around the branch point 62 , as a 
function of z [MeV]. The intersection of the real parts of the two sheets is 
visible. The position of this intersection defines the cut plotted in Fig. 1101 



sociated with the branch point b'^ and indicated with the dotted 
lines in Fig. [12] 

Technically, the transition to the sheets 3 and 4 is achieved by 
following such paths, while further deforming the integration 
contour of the k integration. This further deformation is dictated 
by two requirements; the path in the k plane must not cross the 
poles induced by the denominator of Eq. (fTTT i. and, second, 
must not cross the cut along direction (b) from Fig. |2]in the Zo- 
plane. 

As an example, we show in Fig. [T3] a path in the k plane, 
together with its image in the Za plane given by Za-iz, k) from 
Eq. (fTTT i. With this path, one obtains the continuation to the 
third sheet along the slice s\ from Fig. (TO] shown as the thin 
solid line in Fig. [12] As can be seen in Fig, [13] in the k plane the 

path includes the quasi-two-body singularity at 70 - 400 i 

MeV, while in the Za plane it does not cross the cut of 
indicated with the dashed vertical line. Thus, the integrand is 
always analytic. 

It is instructive to discuss the limiting case of a narrow cr. 
We define this limit by decreasing the crnn coupling constant 
which appears in v""™ of Eq. ©. Then, the mass of the cr 
stays approximately constant. As the cr becomes narrower, the 
branch points b2 and b'^ move towards the real axis according 
to Eq. (flSl l. Simultaneously, the discontinuity on the real axis, 
associated with bi, becomes weaker In the limit of zero width, 
b2 and ^2 coincide on the real axis at z = nit^ + mj]., and the two 
associated cuts coincide as well and run along the real axis. The 
branch point b\ and its cut vanish, i.e., sheet 1 and 2 coincide 
below z - + niff. 

Also, the other sheets overlay: We consider the real parts of 
Ga-N in Fig. [12] They form approximately two "x" with the 
centers at ^2 and bi,. As the cr width approaches zero, the "x" 
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Figure 12: The analytic structure of all four Riemann sheets [arb. units], labeled 
1 to 4, shown close to Re z of the branch points 62 and b'^ . The thick dashed 
lines and the thin solid lines show the structures along the slice si from Fig. 
1101 i.e. slightly below ^2- The thick black solid lines and the thin dashed lines 
are along S2, i.e., slightly above b2. The dotted lines indicate the presence of a 
complex conjugate structure around the third branch point b'^ . 
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Figure 13: Example of a more complicated integration path of the k integration 
in the k [MeV] and the Za- [MeV] plane. The cut in the Za- plane is in direction 
(b), same as in Fig. The result in this example gives the value of Ga-N at 
z = 1600 - 200 (' MeV on the third sheet. 



become symmetric and for zero width overlap exactly, reducing 
the number of Riemann sheets from four to two; the real part 
has then the same structure as in the second row, left, of Fig. [3] 
Thus, in the limit of zero cr width, we have precisely the sheet 
structure of the propagator of stable particles as shown in Figs. 
|2]and|3] with one branch point at z = m'^ + niff, one righthand 
cut, and two sheets. 

In Fig. [T2]we show the labeling of the Riemann sheets. The 
distinction of sheet 1 and 2 is clear: While the physical sheet 

1 is obtained with an integration along an undeformed path in 
Eq. (fTTT i. sheet 2 is the analytic continuation of sheet 1 along 
the righthand ttttN cut. The distinction between sheet 2 and 3 
is defined as follows. We have just argued that a natural choice 
of the cut associated with b2 is along the intersection of the real 
parts (cf. Fig. [TTT i. This also helps us understand the limiting 
case of a narrow cr. Yet, other choices are possible. For prac- 
tical reasons, we prefer a cut that delivers analytic slices along 
the Imz direction, such as plotted in Fig. |9] For this, we have 
to put the cut of i>2 into the negative Imz direction. This is the 
definition which we will adopt for sheet 2: 

Sheet 2, in the lower half plane, is the sheet that can be 
reached along straight paths into the negative Imz direction, 
starting from the real z axis at the continuation of sheet 1 . The 
additional Riemann surface induced by b2, not reachable by 
such paths, is then sheet 3. Sheet 2 in the upper half plane, 
together with the branch point b'^, can be obtained from sheet 

2 in the lower half plane using Eq. (fT3] l. Accordingly, the cut 
associated with b'2 is then in the positive Imz direction. Sheet 4, 
associated with the branch point b'2, is defined accordingly and 
can be obtained from sheet 3 using Eq. (fT3T l. 



Table 1 : Parameterization of the paths in the i-plane for sheets 7 = 2, 3, accord- 
ing to Eq. U6) . The list shows the positions [MeV] of intermediate points in 
the k plane. For the other points kg, k„.-\ , k„^ , see Eq. U7t . Sheet 4 is obtained 
from sheet 3 through Eq. U3t . 
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2.4. Formalism of path deformation 

With these definitions, the paths which lead to the various 
sheets are formally defined. For simplicity, all paths are chosen 
piecewise linear in the k plane. We write for sheet (7); 

G^^(z)= r dklcF{k)g^^^{z,k) 
Jvn 

n2 JJ) , 



2/ 

' = 1 n 



dt (kl - ki^i) [kiit)r F{ki(t)) g'^'iz, ki(t)) 



(16) 



with gaN from Eq. ( fTTT i. Here, \-8a-N^^ evaluated with 
from Eq. Q [n*f^ from Eq. O]. In Eq. (O, A:,(f) = 
ki-\ + {kj - ki-i) t. The k, are the edge points of the paths and 
shown in Tab. [T]for each Riemann sheet { j). For the physical 
sheet j - \,kQ - Q and k„.=i — > 00, i.e. the integration is along 
a straight path from zero to infinity in the k plane. For the sheets 
j - 2, 3, the first, before last, and last point in the k plane are 
given by 



0, 



Ze«{k„.-l,Z) 



2m;, + 100 /MeV forcrAT 
2m;, + 200 /MeV foxpN 
m„ + m/v + 100 / MeV for nis. 



(17) 



where Zeff is the energy of the unstable particle cr,p, or A, given 
by Eq. ( fTTT i for the cr case. For the corresponding expressions 
gpN, gnh as they appear in Eq. ( fT6b for the pN and ;r A propa ga- 
tor, and for the transformations Zp and za, see Refs. liaMJ. 

The distinctions Re z §5 zt, in Tab. [T]ensure that the cut of 
bi is along the negative Im z direction, as discussed at the end of 
the previous section. As an example, we show in Fig. [T4l with 
the solid line the path for the second sheet of the crN propagator, 
case Rez < Rez*, in the Zu plane. Comparing this figure with 
Fig. |7] an additional structure induced by Zo-(^i = 100 MeV) is 
visible; this additional edge point is necessary to avoid mixing 
of sheet 2 and 3. For the endpoint Za-ik^) that is finite in Fig. 
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Figure 14: Integration contour to calculate SG. The plot shows the image 
Za-(k, z) in the Za- plane. See text for further explanations. 



[T4l rather than Zai^i — > as Eq. (fTTT i prescribes, see the 
discussion on implementation in the Jiilich model in Sec. 12.51 



The paths for the sheets j - 2,3 lead to the correct sheets, 
but only in a certain range for the total energy z- Sheet (2), with 
the above defined paths, is calculated correctly for 2m„ + < 
Re z < 2500 MeV and -400 < Imz < MeV, for crN, pN, and 
ttA. For the sheets (3), the continuations are valid around the 
branch points b2', for crN, the range of applicability has been 
tested to be 1600 < Rez < 2100 MeV and -400 < Imz < 
+300 MeV. For pN, 1500 < Rez < 2000 MeV and -400 < 
Imz < MeV; for jtA, 1250 < Rez < 1450 MeV and -250 < 
Imz < MeV. If one wishes to obtain the sheets (3) outside 
these areas, one may have to redefine the paths; the pole in the 
k plane moves as a function of z and for z outside these areas, it 
may cross the integration contour. 

The pN and ttA propagators have the same analytic properties 
as the crN propagator In particular, the poles of the p and the 
A induce additional branch points b2, b'^ of GpN and G„^ in the 
complex z plane. For these channels, these branch points are 
even more relevant for the analytic continuation because the p 
and the A are narrower than the cr and the branch points lie 
close to the real axis according to Eq. ( fTSI ). The branch points 
b2 and b'2 have large numerical effects in their surroundings as 
can be seen in Figs. |9] to [12] and, thus, play an important role 
for the analytic continuation. The integration paths to obtain 
the different sheets of the pN and nh. propagators are given in 
Table[Tl together with the points defined in Eq. ( fTTl i. 

To finish this discussion, let us again point out the importance 
to control the cuts induced by the branch points b2, b^. As dis- 
cussed above, the condition Re z §i Re Zhi controls the cut. In 
particular, in these two different cases the k integration contour 
passes by the quasi-two-particle singularity above and below in 
the k plane, respectively. This distinction is mandatory; if the 
contour passes for all z above or for all z below that singular- 
ity one evaluates the third sheet instead of the second for some 
z; poles on that sheet are much less relevant than poles on the 
second sheet, because the third sheet is not directly connected 
to the physical axis. 



2.5. Implementation of the continued crN, pN, ttA sheets 

To obtain the amplitude of the Jiilich model on the different 
crN, pN, ttA sheets, we add the difference between second and 
first sheet, similar as in the case of the channels ttN and rjN in 
Eq. ([8]i. In Sec. 12.11 we have seen that for the ttN, tjN prop- 
agators the prescription to obtain the second sheet consists in 
adding the discontinuity twice. In the present case of the un- 
stable crN, pN, nA propagators, one can proceed similarly. Yet, 
for the effective ttttN channels, this implies an approximation 
which will be discussed below. 

In Fig. [14] the image Zo-(^, z) of the A:-integration in the Zo- 
plane is shown. Consider the case j - 2 for the crN propagator, 
i.e. the second sheet. The dashed line shows the integration 
contour for j - I, the first sheet. In the k plane, this integration 
is along a straight path from = to A: = 00. 

The soUd lines show the integration according to the case 
7 = 2 with the edge points ko, ki, and k2 given in Eq. ( fTTI ) 
and Table [1] For the endpoint k^, we have chosen, instead of 
kj, — > 00, an intermediate point Za-ikj.) that is on the contour 
of the integration for the first sheet and below the nn threshold 
(remember that for the second sheet, the nn cut is in direction 
(b) of Fig. fT4l) . In the difference of second and first sheet, the 
remaining path from the point k^ to 00 cancels. 

For the plots of numerical results in this section, we intro- 
duced a form factor F in Eq. (fTTI) to regularize the integral. 
For the implementation in the Jiilich model, this artificial form 
factor is removed; the result for the difference of sheets is finite 
anyways as discussed before. 

Thus, the difference between sheet (j) and sheet (1) is given 

by 



Geif)f=l, 



(18) 



where j - 2, 3. Here and in the following, the subscript "eff" 
indicates the effective nnN channels. The end points of both 
integrations for G^J^ and G*|j! are at 

Zeff(k„.,z) 

_ f 2m;r - 1 00 Me V + ;■ Im z {o-N, pN) 

(m;r + wai - lOOMeV + ;Tmz (nA) 



(19) 



which replaces k„. from Eq. ( fTTI i. As Fig. [14] shows, the 
overall integration path to obtain 6G is a closed contour, 6G - 
j> dk f{k,z), but with start and end point on different sheets of 
the self-energy of the unstable particle. 

The next step to implement the sheets of the effective nnN 
propagators in the Jiilich model is similar to Eq. ([8]l, 



6G^^ 



r , 2 (qcd\V\qeff}{qeff\T'^'Hqah} 

I dq^sq^ff 



Z - £eff - Ileff 



^G^^ = S&il{qcd\V\qU^«)){qlJ.^nT'%ab) 



■Mj) 



(20) 



As in case of the channels with stable particles, 6G is added at 
the on-shell point. In case of the channels with stable particles, 
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this is exact because the 6 function that evaluates the imaginary 
part puts the vertices automatically on-shell [cf. Eq. (|7]i] . 

In case of the effective nnN channels, this factorization of 
the last line in Eq. ( |20] | is an approximation. For q^^{tW), we 
choose the on-shell momenta of the kinematics 

^^„(eff) = qlJ^z -^rriN+n m^) (crN, pN) 

q^nieS) = ^on(z ^m„ + niA) (ttA) (21) 

with mA = 1232 MeV and m^, m„ the nucleon and pion mass, 
respectively. We test the cases « = 2, 3 for all results. The pole 
positions and residues are very similar in both cases. 

Another motivation for the factorization is given by the fol- 
lowing considerations. The kinematic z + 2 m„ for the 
(xN (and also pN) propagator corresponds in good approxi- 
mation to the maximum of Imk^ga-N from Eq. ( fTTT i. which 
gives the discontinuity of the nnN cut along the physical axis. 
Strictly speaking, the discontinuity receives contributions from 
Imk^ ga-N for all k, but the distribution is still concentrated 
around the maximum due to the factors k- and the form factors, 
that suppress the contributions for small and large k values, re- 
spectively. 

Furthermore, the kinematics z — > niN +2 m„ or z — > 111^+3 m„ 
corresponds to typical 2n invariant masses in nN — > nnN in the 
second resonance region in nN nnN, where the branching 
ratios into pN and crN are analyzed. 

The principal difficulty in going beyond the factorization of 
5G*^ in Eq. ( |20] | is the incompatibility of the integration paths 
that lead to unphysical sheets of the effective nnN propagators, 
and the position of three-body cuts in the transition potentials V 
that require a diff'erent integration path. While the access to the 
different sheets is dictated by the deformed paths as discussed 
above, the momentum integration that respects three-body cuts 
is from k - Q along a straight path rotated into the complex 
plane ll20ll . A unified treatment is beyond the scope of this 
work, although possible in principle. 

The analytic continuation of T for the first sheet, as given by 
the integral terms of Eqs. (jS) and ( l20t . is restricted by the two- 
body cuts of the stable nN and riN propagator, as well as the 
pseudo-two-body cuts of crN, pN and nA. These zeros of the 
propagator denominators z- Ei- E2, together with the rotation 
of the integration path into the lower complex q half plane lEoll . 
induce fallacious non-analyticities in the lower z half plane. 
Thus, we restrict the analysis of the analytic continuation to 
the upper z half plane. Yet, for all results of this study in Sec. [3] 
we will quote the values for pole positions, residues etc. in the 
lower z half plane that are easily obtained using Eq. ( fT3T l. 

Additionally, the three-body cuts of the potentials V (in par- 
ticular, the pion exchange), together with the rotated integration 
path, induce similar structures for the first sheet in the upper z 
half plane. Since these structures are located at Imz > 150, 200 
MeV for all z we can safely search for resonances up to a width 
of r = 300 or 400 MeV, and this issue is of no relevance here. 
For energies above the second resonance region, even larger 
values for Imz are accessible. Within our formalism, an ex- 
change propagator is of the form 1 /[Ey-{z - E^ — E\ — E2 + ie)}. 
There are kinematic factors IjE^ of the exchanged particle x. 



They induce similar structures as the three-body cuts but they 
are also situated at Imz > 150, 200 MeV. We make no attempt 
to access the sheets induced by the three-body cuts associated 
with the exchange potentials V. These sheets are far away from 
the physical axis and structures on those sheets can only have 
minor impact on physical scattering. For completeness, the 
short nucleon, circular and other cuts, associated with nucleon 
exchange and crossing symmetry, are discussed in AppendixlAl 
It is shown that this cut structure is indeed present in the model 
and can be identified with different ingredients of the model. 

3. Results 

3.1. Selection of Rietnann sheets 

In Sec. |2]the analytic structures of the propagators of chan- 
nels with stable particles nN, rjN and of effective nnN channels 
ttA, pN, and crN have been determined. The analytic contin- 
uations of the propagators determine the analytic continuation 
of the T matrix. In this section, we determine the properties of 
T{z) in terms of poles and zeros in the complex z plane. 

For a channel with stable particles, there are two sheets, 
while for unstable particles, there are four as we have seen in the 
previous section. Thus, for the channel space considered here, 
there are 2^4^ = 256 sheets corresponding to the two stable and 
three effective nnN channels. 

Poles can be located on all sheets. Yet, depending on the 
sheet, their influence at the real z axis ("physical axis") is differ- 
ent: For a channel of stable particles with two sheets, consider 
a virtual state. This is a pole at z = zo on the second ("unphys- 
ical") sheet in the lower z half plane, and below threshold, i.e. 
Rezo < Zthres- Then, the closest point to the physical axis is 
Zthres [cf. Fig. [3] third row]. Because the pole contribution has 
a 1 /(z - Zo) energy dependence, such a pole can only appear as 
a structure at Zthres, e.g. as a threshold enhancement. It cannot 
create a resonance shape on the physical axis at z = Re zo- 

Thus, for the pole search with respect to the nN and rjN chan- 
nels, below threshold one chooses the first sheet of the propa- 
gator G. This means looking for bound states with respect to 
that channel. Above threshold, one chooses the second sheet. 
Formally, this can be written as 

^^|g(2) ifRez>Zthies 
|g(i> ifRez<Zthies 

(22) 

where Zthies - m + M is the threshold energy. This choice 
will induce cuts in the amplitude along the Imzo direction, for 
z — m„ + and Rez — + m^. 

As pointed out in Sec. 12.31 for each effective nnN channel 
nA, pN, and o-N there are four sheets. For the pole search, we 
choose the unphysical sheet with respect to the nnN cut, i.e., the 
second sheet according to the definitions at the end of Sec. 12.31 
In particular, the second sheet is reachable from the physical 
axis via straight paths in the negative Imz direction. There are 
additional sheets induced by branch points b2, b'^ in the com- 
plex plane as we have seen in Sec. 12.31 Yet, coming from the 
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Figure 1 5 : Modulus | r(S 1 1 )| [arbitrary units] as a function of the complex scat- 
tering energy z. [MeV]. The poles of the N'(i535) and A'*(1650) are clearly 
visible. Also, one sees a discontinuity along Rez = m^, + m^, coming from the 
prescription of Eq. )22t for the rjN channel. Behind the ^'(1650), one sees a cut 
lying in the positive Imz direction, induced by the branch point zi/^ = + Zp 
of the pN channel. 



compared to the a terms, and the replacement G^^ — > G^^y does 
not change much the position of the zero; the resonance pole 
will reappear on the third or even fourth pN sheet. 

Such replica of poles on other sheets have no physical impli- 
cations. E.g., the A*(1700) has a pole atzo = 1637 - 118 ;MeV 
on the second pN sheet and a replica on the third one, just a few 
MeV away from zo and with almost the same residues. Even 
the distance from both pole positions to the physical axis via 
paths over analytic regions of the amplitude is approximately 
the same. Yet, one of the poles is sufficient to describe the A* 
properties on the physical axis to high accuracy liSftl . The two- 
pole structure is, in this case, trivial. 

The situation is different for e.g. the proposed two-pole struc- 
ture of the A(1405) ll33ll . In the latter case, the two poles lie 
on the same sheet on different positions and have quite differ- 
ent residues to the different channels. While there is evidence 
[SSA for the two-pole structure of the A( 1405), the trivial replica 
found here have no physical consequences. 

Yet, in the present context, we sometimes find a pole on the 
third sheet of pN, but no counterpart on the second sheet; in 
that case the coupling of the state to pN is strong. In such a 
situation, discussed in Sec. 13.61 a pole on the third pN sheet 
can indeed induce visible structures on the physical axis that 
cannot be explained from the amplitude on the second sheet. 



physical axis, those sheets are only accessible via paths around 
these branch points and poles on them have little influence on 
the physical axis [cf. Fig. (TO). This is in analogy to the case of 
the ttN and rjN propagators discussed previously. 

The selection criteria discussed above define one out of 256 
sheets, where poles are searched. We refer to this as "second 
sheet" of T in the following. As an example of the analytic 
structure. Fig. [13] shows |r| on that sheet, in the S n partial 
wave in ttN — > ttN. The two poles associated with the A^*(1535) 
and the A^*(1650) are clearly visible. Yet, there are additional 
structures: First, a discontinuity along Re z = m,j + coming 
from the prescription of Eq. ( |22] |. Second, one sees a similar 
discontinuity starting behind the A^*(1650) pole resulting from 
the definition of the second pN sheet: the branch cut is in the 
positive Imz direction and starts at the branch point ziy^ - + 
mn ~ 1700 + 64 ;■ MeV where Zp - 763 - 64 / MeV is the pole 
position of the p [cf. Eq. (flSTll. 

The first sheet is free of poles as we have checked; we search 
for poles on the second sheet. Yet, there are poles on other 
sheets which we comment on in the following. Consider a pole 
that couples only weakly to a given channel, i.e. its residue 
to e.g. the pN channel is small. Suppose the pole has been 
found on the second sheet. Then, at the pole position, the term 
(1 - VG) ' which appears in the solution of Eq. (H} is singular. 
We consider an element of this matrix and write symbolically, 
with G^^ the pN propagator on the second sheet. 



(l-VG)^a + bGf^^O 



(23) 



omitting further indices, sums and integrations, a contains the 
terms with intermediate states of other channels. The weak cou- 
pling to pN is reflected by the fact that the b terms are small 



3.2. Pole positions and residues 

The scattering amplitude t for the transition / ^ / in channel 
space is connected to the amphtude T of Eq. ([U by 



kEo) 



(24) 



where k {E, cj) are the on-shell three momentum (baryon, me- 
son energies) of the initial or final meson-baryon system. In 
order to extract the pole residue, we expand the amplitude T*^' 
on the second sheet in a Laurent series around the pole position. 



-I 



z-zq 



+ ' +0(_z- zo)- 



(25) 



The residue a_i and constant term ao can be obtained by a 
closed contour integration along a path r(z) around the pole 
position Zo, 



an = ^.(£ 



P^\z)dz 
r(z) (z - zo)"+' ■ 



(26) 



Alternatively, a_i can be expressed in terms of dressed quanti- 
ties H, 



-(t) 

D 



1 - #2 

az 



(27) 



where Yd (r^\ ^) is the dressed annihilation vertex (creation 
vertex, self-energy) as defined in Ref. ifsoll . evaluated on the 
second sheet at zo- We have explicitly checked for all reso- 
nances, using the explicit values for To, 2, that the results of 
Eqs. ( |26] | and dZTb agree. 
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The dressed quantities in Eq. ( |27] | are given in terms of the 
non-pole part T^^, which dresses the bare creation and anni- 
hilation vertices jb- For details see Ref. ifsoll . Note that 
for a simple energy and momentum independent s wave inter- 
action, y^" - jB while for higher spins and partial waves the 
connection between bare annihilation and creation vertices can 
be more complicated. 

Using Eqs. (l24l i and (IZST i, the pole residues R - \R\e'^ as 
quoted by the PDG |01 can be calculated. For the residue phase 
we consider the usual ifsill definition given by 



|/?|< 



M-z- ir/2 



(28) 



for a resonance with width F on top of a background tb- Com- 
paring Eq. ( |28] | with Eq. (IZSl l and using Eq. (l24l i. the pole 
residue R and its phase are given by 



\R\ = \a-ip„ 



- -n + arctan 



Im (fl-lP;riv) 



Re {a-ip„N) 



(29) 



where p„N is the phase space factor p from Eq. (l24l) for the 
ttA^ — > ttA^ transition, evaluated at the complex pole position. 

Poles of the amplitude are searched for on the second sheet, 
as defined and described in Sec. |2] The results for pole 
positions and residues are summarized in Table |2] The ex- 
tracted resonance parameters are compared with other studies 
mH,!!, all of them accepted by the PDG M- 

For prominent resonances with large branching to ttN, the 
different analyses are in reasonable agreement with the present 
results. For other resonances that are wide and/or couple only 
weakly to nN, the results are much more disperse and there are 
noticeable differences among the results of Refs. fs, 5^ ^ and 
also to the results of the present study. Note that for resonances 
such as A*(1910) and A^*(1720), little is known about residues 
and phases. Given, e.g., the span of 6 - -1-172'', -90° from 
the PDG for the A*(1910), it is no surprise that the value of 
the present study of 6 - -153" does not match any of the two 
results. 

The Roper A^*(1440)Pii resonance does not require a gen- 
uine pole term in the Jiilich model jl9]. Instead, the resonance 
shape is dynamically generated from the coupled channel inter- 
action together with the unitarization from Eq. ([T]). Here, we 
can confirm this result, because we have indeed found a pole 
for this resonance [cf. also Sec. 13.311 . 

A special situation is given for the S 1 1 partial wave in which 
two resonance interfere making the extraction of resonance pa- 
rameters more difficult. Table |2] shows that the values from the 
PDB for both \R\ and are very different in the various stud- 
ies for the A^*(1535). This is due to the systematic uncertainties 
from the close-by rjN threshold plus the interference with the 
A^*(1650). The values of the present study lie within these wide 
spans. This issue is further discussed in Sec. 13.51 

For the other resonances, the results of the present study 
sometimes lie at the borders of the ranges quoted in the PDG. 
For the A(1232), the present values for pole position and residue 
lie even outside the range of the other studies. In Fig. [16] the 



Table 2: Resonance pai'ameters in the present study. The zo are the pole posi- 
tions. The moduli and phases 6 of the residues correspond to the nN decay 
channel. For every resonance, the first line quotes the results of the present 
study, followed by the values from Refs. Sp^lfijl as quoted in the PDB Q|. 
Resonances are included in the Jiilich model via explicit .v channel exchanges 
except for the Roper resonance which appears dynamically generated. 
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Figure 16: Imaginary part of tlie P33 amplitude witli tlie A(1232). Tliere are 
some residual deviations of the present fit (red solid line) and the SES [FA08] 
from Refs. (57, J8J (data points). 



Table 3: Position of zeros of the full amplitude T in [MeV]. There is always 
another zero at the complex conjugate position according to Eq. )13t . For 
comparison, the zeros determined in Ref. (sgll (FA02) ai'e also quoted. 



first sheet 
Pii 1235- 
D33 1396- 



second sheet 

0/ Su 1587-45/ 

78/ 531 1585 - 17/ 

P31 1848 - 83/ 

Pn 1607-38/ 

^33 1702- 64/ 

£)i3 1702-64/ 



Ref. [59j 
1578 - 38 / 
1580-36/ 
1826- 197/ 
1585 -51 / 

1759-64/ 



present solution is shown together with the SES solution of the 
partial wave analysis of Ref. 15711 . There are residual discrep- 
ancies for Im P33 which suggest that the Jiilich model may need 
some fine-tuning to obtain an improved fit. It is a very small 
effect, but the PWA results are quite precise for the A(1232). 
Note, however, that even for near identical phase shifts K ma- 
trix analyses and analyses using analytic T matrices will give 
in general diff'erent pole positions and residues as a result of the 
different analytic properties. 

Poles and zeros are important parameters and determine the 
global appearance of a partial wave amplitude. In the one- 
channel case, unitarity leads to zeros on the first sheet if there 
is a pole on the second sheet. For the coupled channel case dis- 
cussed here, there is no such direct connection between poles 
and zeros. Yet, zeros play an important role. For example, as 
pointed out in Ref. |20(]. the unitarity constraint in nN to rjN 
leads to a zero above the A^*(1650) in the absence of additional 
inelasticities; this results in an unobserved dip in the nN t]N 
cross section. The introduction of couplings of the S \ \ reso- 
nances to ttA in Ref. lEoll solved the problem, simultaneously 
reducing the maximum of the nN rjN cross section at the 
A^*(1535) energies to the physical value. 

The zeros of the Jiilich model have been determined in the 
present study, with the results shown in Table |3] There is a 
zero on the physical axis for the Roper channel. This reflects 



the peculiar phase shift of the P\i partial wave as discussed in 
Ref. t2ft1 . For the 531, P31, and f 13 partial waves, the various 
zeros on the second sheet are in close vicinity to the respective 
resonance poles. The zero of 5 n lies in between the two Sn 
resonances and will be further commented on in Sec. 13.51 

In Table [3] also the zeros extracted from the FA02 solution 
of Ref. i59r are shown. In their sheet numbering, our sheet 
2 is their sheet 1. The global pattern is similar; in S'n the 
zero is in between the two resonances, in 531, P31, and P13 
the zeros are correlated with the respective resonance poles. 
Most interestingly, we find zeros in P33 and D13 at the pN 
branch point associated with the quasi-two-particle singularity 
at ZbXpN) = 1702 - 64 / MeV. 

3.3. Pole structure of the Roper resonance 

The poles from Table |2] all lie on the second sheet. We have 
also searched for poles on other sheets for some selected cases. 
E.g., in the partial wave analysis SP06 of Ref. [5], a pole of 
the Roper has been found at zq - 1359 - 81 / MeV on the sec- 
ond ttA sheet (in their counting: first sheet), and another one 
at zo = 1388 - 83 / MeV on the third ttA sheet (their count- 
ing: second). Also in the present study, we find a second pole 
of the Roper on the third ttA sheet at zq - 1387 - 71 / MeV 
which is just a few MeV away from the pole on the second 
sheet quoted in Table |2] As discussed in Sec. I3.1l this is rather 
a replica of the pole on the second sheet, without physical im- 
plications, than a genuinely new structure; indeed, within the 
Jiilich model the coupling strength of the Roper to the nA chan- 
nel is moderate ||60ll ; a change of sheets does not change much 
the resonance properties in this special case. This is also re- 
flected in the value of the residue of the additional Roper pole: 
a_ I (2nd 7rA)/fl_i (3rd TT A) = 1.06 - 0.01 / for the ttA^ residues of 
the poles on the second and third nlS. sheets. 

The rather different pole positions and residues of the two 
Roper poles in Ref. 151] indicate a larger coupling to the ttA 
channel, and in this case the branch point ^2 plays an important 
role as stressed in Ref. jsl which makes a simple Breit-Wigner 
parameterization of the Roper questionable ISfl. 

In the Jiilich model, the nnN inelasticity of the Roper is 



rather given by the effective crN channel 0181 11911 . The early 
onset of inelasticity in the P\i partial wave is naturally ex- 
plained by the s wave character of the crN coupling to the 
Roper ifisl [19I1 . In contrast, tt A couples in p wave to the Roper, 
and the centrifugal barrier renders the contribution small at low 
energies. Thus, a large ttA coupling to the Roper would be 
needed to provide the necessar y in elasticity at low energies. 
Also, it has been shown in Ref. M19ll . that the persistently large 
inelasticity at higher energies can be explained more naturally 
by a large crN coupling. 

For completeness, let us mention a similar situation for the 
pN sheet: on the third pN sheet, there are poles in P13 and D33 
situated a few MeV away from their counterparts on the second 
pN sheet, associated with the A^*(1720) and A*(1700), respec- 
tively. Again, none of the two resonances couples strongly to 
pN 16011 . and replicas of the poles on the second sheet appear 
on the third sheet. The pN sheet in D13, however, is special and 
will be discussed separately in Sec. 13.61 
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Table 4: Resonance couplings gj [10 ^ MeV to the channels nN and rjN. 





nN 




j]N 


N*{l440)Pu 


11.2- 


5.0/ 


-0.1+0.0/ 


N*(1520)Di3 


8.4- 


0.8/ 


0.16-0.60/ 


A^*(1535)5ii 


8.1 + 


0.5/ 


11.9-2.3/ 


A^*(1650)5ii 


8.6- 


2.8/ 


-3.0 + 0.5/ 


A^*(1720)Pi3 


3.7- 


2.6/ 


-7.7 + 5.5/ 


A(1232)P33 


17.9- 


3.2/ 




A*(1620)53i 


2.9- 


3.7/ 




A* (1700)^33 


4.9- 


1.0/ 




A*(1910)P3i 


1.2- 


3.5/ 





Table 5: Bare and renormahzed vertices y'- and V'- in [10^^^ MeV"''^] for some 
resonances in the nN — > nN transition. The last two columns show the ratios 
defined in Eq. ^31^. 



and dressed vertices, as defined in Ref. ifsoll . The other columns 
show the ratios 

r = KTd - 7B)/rD\, 

r' = |1 - Vl -S'l, (31) 

i.e. the relative differences of dressed vertex versus bare ver- 
tex, and coupling versus dressed vertex. The renormalization 
of the bare vertex (ratio r) is large in almost all cases. But also 
dressed vertices and couplings can be very different (ratio r'), 
due to the appearance of the term 1 - Z' in the denominator 
of the expression for the residue from Eq. dZTl i. The conclu- 
sions from this behavior are the same as in Ref. llsoll : bare 
and dressed vertices are model dependent quantities; the non- 
pole T matrix r"^^ from the model dependent decomposition 
T = + T^^ enters in the calculation of To- In contrast, the 
couplings gi provide a meaningful expansion parameter of the 
amplitude around the pole of the resonance, independent of the 
amplitude decomposition into pole and non-pole parts. 









r[%] 


r' [%] 


A^*(1520)Di3 


6.4 - 0.6/ 


13.2+1.2/ 


53 


61 


N*{n20)Pi3 


-0.1+5.4/ 


0.9+4.8/ 


24 


45 


A( 1232)^33 


1.3 + 13.0/ 


-2.8 + 22.2/ 


45 


40 


A*(1620)53i 


0.1 + 14.3/ 


5.0 + 5.7/ 


130 


66 


A*(1700)D33 


5.4 - 0.8/ 


6.7+ 1.0/ 


33 


54 


A*(1910)P3i 


9.4 + 0.3/ 


1.9-3.2/ 


222 


22 



3.4. Couplings 

It is convenient to express the residues a_i in terms of a few 
parameters g, given the different residues for the transitions 
within n channels. It is possible to write, for the residues into 
the 7tN and rjN channels. 



(30) 



with a unique set of g, quoted below. The g, will be referred to 
as couplings in the following. Note that in the determination of 
the gi there is an overall undetermined sign which we have fixed 
by choosing the real part of the coupling constant g„N positive. 

In Table |4] we list the coupling strength g, of the resonances 
to the nN and rjN channels. The couplings to the effective nnN 
channels will be published elsewhere i60ll . They are important 
in the description of the nN — » nnN reactions. 

The allowed couplings gi in Table |4] are all non-zero, while 
only a few bare couplings are included E.g., the bare tjN 
coupling of the A^*(1650) is zero while the corresponding g in 
Table |4] is finite. This is because the rescattering in the unitary 
coupled channel model renders the residue finite even when the 
bare couplings may be zero. 

In Ref. [SOD, the bare and dressed vertices and the residue 
of the A(1232) have been evaluated. In Table |5] we quote the 
corresponding results for other resonances. The vertices have 
been evaluated on the second sheet at the pole positions of the 
respective resonances. The second and third columns show bare 



3.5. Resonance interference in S n 

As an example of the analytic structure we discuss the S \ i 
partial wave. This is of particular interest, because there are two 
resonances in this partial wave. For a systematic discussion of 
other partial waves and their analytic structure, in particular Pn 
with the Roper resonance, see Ref. ll60ll . 

In Fig. [17] the amplitude r for the S \ \ partial wave is plotted, 
obtained from the amplitude T via Eq. (|24] |. The full solution 
of the Jiilich model is indicated with the solid red lines. It de- 
scribes well the SES solution of Ref. 115 711 up to z ~ 1 .9 Ge V. 
The blue dotted lines indicate T^^ . We can also plot the pole 
approximation from Eq. (l25T l. For simplicity, we set qq - 0. 
On the physical axis, the pole approximations of the A'*(1535) 
and A^*(1650) appear as resonant like structures indicated with 
the black dashed-dotted lines. We first discuss the amplitude 
above the rjN threshold. For a discussion of the cusp and the 
amplitude below the rjN threshold, see below. 

At first sight, the shapes of the two resonances in Re t are 
quite different: While the pole approximation of the A^*(1535) 
shows a familiar shape with a maximum and a minimum in 
Re 5 11, the A^*(1650) looks quite different. The reason is that 
a_i is a complex number that mixes real and imaginary parts 
of a classical Breit-Wigner shape. In other words, the phase of 
the resonances from Table |2] is responsible for this twisting of 
resonance shapes and can have a very large effect. 

The individual contributions from the two resonances (black 
dashed-dotted lines) are quite different from the full solution. 
However, the sum 



a 



z-z, 



-1 

^•(1535) 
■0 



a 



Z-Z, 



-1 

A" (1650) ' 




(32) 



indicated as the red dashed lines, fits the full solution quite well 
over the entire resonance region. Thus, the two resonances 
cannot be treated separately but must be treated together; the 
residue from the A^*(1535) provides a strongly energy depen- 
dent background in the A^*(1650) region and viceversa. Note, 
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Figure 17: Amplitude in the 5ii partial wave. The data points represent the 
single energy solution from the partial wave analysis [FA08] of Refs. I'ST'.'sI. 
The solid red and dotted blue lines show the full amplitude and T^^, respec- 
tively. The dashed dotted lines show the pole approximations for the A'*(1535) 
and the A'*(1650) (from physical poles above the t]N threshold, from hidden 
poles below). The red dashed lines show their sum. 



e.g. for Rer the strong energy dependence of the tail of the 
A^*(1535) in the region of the A^*(1650). The resonances inter- 
fere which each other. 

Thus, for a theoretical description one needs a unitary cou- 
pled channel model like the present one, which also allows for 
resonance interference. Otherwise, if one tries to extract res- 
onance parameters individually for each resonance, one needs 
a substantial phenomenological background. Then, the param- 
eters depend very strongly on that particular background and 
results are not reliable. 

At the level of pole positions and residues, resonance in- 
terference appears as a strong cancellation effect of pole ap- 
proximations as discussed before; this should be clearly dis- 
tinguished from what is conventionally meant with resonance 
interference: in a microscopic approach like the present one, 
the full partial wave amplitude with more than one explicit res- 
onance can be decomposed according to 

T — T^^ + 

- Z V ^^o)r' (33) 

with resonance indices r, r', dressed vertices To and bare prop- 
agator 5 b as defined in Ref. for the one-resonance 
case. Resonance interference is allowed by non-vanishing off- 
diagonal self energies S12 and 221. Second, even for S12 = 
2^21 - 0, the individual parts of the sum in Eq. (1331 ) can be 
large and make it difficult to phenomenologically disentangle 
the resonances, as in case of the A^*(1535) and A^*(1650) dis- 
cussed here. 

Note that the full amplitude T in Eq. ( l33b is unitary which 
is automatically ensured by the complex phases of the dressed 
vertices F, that are linked to T^^ [cf. Ref. ifs^ l: the expres- 
sion in Eq. ( l33b for should not be confused with a sum of 
two Breit-Wigner resonance amplitudes. While a single Breit- 
Wigner amplitude is unitary, a sum of two is not. 

In contrast to the decomposition of Eq. ( |33] ), the expansion 
of the amplitude in terms of poles and residues, as given in 
Eq. ( l32b . does not rely on the model dependent separation into 
and T^^ ifsoll : Eq. ( l32b allows to study the interference of 
resonances independent of the decomposition. 

For the pole search it is important to find all poles, even if 
they couple only weakly to ttN. For example, in Ref. ifsoll we 
have found poles in the so-called non-pole part of the amplitude 
T^^ (in P33) which are dynamically generated and have no gen- 
uine pole term explicitly appearing in the amplitude. A useful 
way to identify poles and zeros is to consider the two curves 
defined by Re r*-'(z) = and Im T^^\z) - 0. In such a "Gauss 
plot", the curves intersect at poles and zeros. 

In Fig. [18] the analytic structure of the amplitude r^^'(z) 
on the second sheet is visualized in a Gauss plot. The solid line 
is defined by Re T^^\z) - 0, the dashed line by Im T^-\z) - 
0. The two poles from Table |2] and the zero from Table [3] are 
indicated in the figure and lie indeed at the intersections. There 
are no further intersections in the resonance regions. 

The dashed dotted and dotted lines show the Gauss plot ob- 
tained from Eq. (|32] | instead of the full amplitude. The sum 
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Figure 18: TheSii amplitude T'-^\z) on the second sheet. The lines are defined 
by Re T'-^H.z) = and Im T'-^Hz) = 0. See text for further explanations. 



from Eq. ( l32l l reproduces remarkably well the amplitude even 
further away from the pole positions, including the position of 
the zero in between the two resonances. The kinks of the full 
solution (dashed line) at Re z = niij + nif^ and Re z = Re zt', ipN) 
originate from the cuts appearing in the definition of the second 
sheet as discussed in Sec. 13.11 

This comparison of full solution and the sum of Eq. ( l32l i. on 
the real axis and in the complex plane, reflects the accuracy of 
the extraction of pole positions and residues carried out in this 
study. There are, of course, still residual deviations between 
the full solution T^^\z) and T^a\z) from Eq. (|32] |. e.g. for the 
position of the zero. They are due to higher order terms in the 
Laurent expansion which are, however, relatively small for the 
two 5 11 resonances. 

In the following, the amplitude below the rjN threshold is 
discussed. This part of the physical axis is directly connected 
to the part of the complex plane where the nN channel is on the 
second sheet (2) but the rjN channel is on the first sheet (1). This 
sheet is called 21 in the following. The above discussed pole 
structure is not on sheet 2 1 , but on sheet 22 (directly connected 
to the physical axis above the rjN threshold). 

Thus, the discussed pole structure on sheet 22 has no im- 
pact to the physical axis below the tjN threshold. However, on 
sheet 21 itself, there are also two poles. The real parts of their 
pole positions are greater than m^, + rriN- These "hidden poles" 
thus lie in the same z region as the discussed A^*(1535) and 
A^*(1650), but on a different sheet. This situation is illustrated 
in Fig. [19] The blue arrows indicate how the poles on differ- 
ent sheets are connected to the different pieces of the physical 
amplitude (red solid Une). The figure shows schematically the 
different poles of the A'^*(1535); for the A^*(1650), not shown in 
the figure, the situation is the same. 

For the hidden poles on sheet 21, one can also draw the pole 
approximations and their sum from Eq. ( l32l i. This is indicated 
in Fig. [17] for m„ + < z < + m^, with the black dashed 



N*(1535), sheet 22 




Re z [MeV] 



Figure 19: Hidden poles in the 5 n amplitude. The A^'CISSS) and A'*(1650) on 
sheet 22 {ttN, rjN second sheet) are responsible for the resonant shapes above 
the r]N threshold, the hidden poles on sheet 2 1 (;rA' second, rjN first sheet) are 
visible through their shoulders on the physical axis below the rjN threshold. 
The arrows indicate the impact of poles on different sheets to different pieces 
of the physical amplitude. 



dotted lines and the red dashed line. 

As figure [17] shows, the pole approximations from above 
the rjN threshold, in combination with the pole approximations 
from below the rjN threshold, provide a good reproduction of 
the pronounced rjN cusp (note T'^^ shows no cusp effect at all). 

It is also worth mentioning that all resonance contributions 
become small as one approaches the nN threshold (cf. dashed 
line in Fig. [TTI ). This is a consequence of chiral symmetry 
which is included in the Jiilich model through the modifications 
of Ref . ll20ll . There, a derivative ttNN* (1650) coupling provides 
the disappearance of resonance contributions close to threshold 
as required by chiral symmetry. 

3.6. A pole on the third pN sheet 

The analytic properties of an amplitude with unstable parti- 
cles can imply complex structures one of which is discussed in 
the following. In Fig. [20] the D13 partial wave is shown. The 
full solution is indicated with the red solid lines and reproduces 
well the partial wave analysis from Ref. 157 1. 

The blue dotted lines represent the non-pole part T^^ as de- 
fined in Ref. |50], i.e. the amplitude without the ^-channel 
resonance exchange diagrams. The pole term contains 
the s channel exchanges and the full amplitude is given by 
T - T^^ + ifsoll . A resonant structure in r"^^ is visible 
at around z ~ 1.7 GeV, which disappears in the full solution 



T - T^^ + T^. The second Riemann sheet of T^^, however, 
is free of poles. Instead, the cut structure from the pN branch 
point b2 at Zbi - 1702 - 64/ MeV from Fig. [TO]is enhanced. 
This is a sign that there is structure on the third pN sheet [cf. 
Fig- [12)1 . This is indeed the case. Using the prescriptions from 
Sec. 12.41 it is possible to analytically continue the amplitude of 
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Figure 20: Amplitude in the D13 partial wave. Labeling of the plot as in Fig. 

[m 

the Jiilich model to the third pN sheet. Indeed, there is a pole 
at 1613 - 83 / MeV. It has a strong coupling to the pN channel 
and a medium size coupling to nN and is, thus, a state dynami- 
cally generated mainly from the attractive interaction in the pN 
channel. 

Fig. l2n illustrates the effects of a pole on the third sheet: A 
path from the pole position to the physical axis must necessarily 
pass around the branch point b2 to get from the third sheet to the 
second sheet that connects to the physical axis. The part of the 
path on the third sheet is indicated with the dotted lines in Fig. 
\2l\ The straight red line indicates our choice of the cut, where 
second and third sheet are connected [cf. Sec. 13.111 . Once on 
the second sheet (solid blue line), it is possible to approach the 
physical axis on a straight path. 

In other words, the pole on the third sheet can affect the phys- 
ical amplitude only via a detour around b2- The structure in the 
physical amplitude is, thus, rather located at the position of b2 
than at the actual pole position. 

Summarizing, the branch points b2 of effective nnN channels 
lead to complex structures of the amplitude. The structure that 
appears as a resonance on the physical axis is in fact a threshold 
effect of the quasi-two-particle threshold pN associated with b2', 
it is induced by a pole at a different position, on a sheet that is 
not directly accessible from the physical axis. 

As a test, one can decrease the pnn coupling of the p self- 
energy. As a consequence, the p becomes narrower and b2 ap- 
proaches the physical axis. Then, the structure that seemed to 
be a resonance in r"^^ for the physical p, appears as a typical 
threshold cusp on the physical axis. 



Figure 21: The D13 pole in T^^ on the third pN sheet and its influence at the 
physical axis. The inset shows schematically the structure visible in T^^. 



The mechanism discussed here may also apply to a very sim- 
ilar structure observed in Ref. 116 ill . There, AK scattering and 
its implication for the pentaquark were discussed. For finite 
A width, a resonant structure appears close to the AK thresh- 
old which turns into a cusp structure once the A width is set to 
zero, in analogy to the present case; thus, what was noticed as a 
fading away of the pole ll6lll for some combinations of input pa- 
rameters, may correspond to the pole moving far into the third 
sheet, in the present formulation. 

Coming back to the discussion of Do in the Jiilich model, 
in the full amplitude ifsoll T — T^ + T'^^ , the dynamically gen- 
erated pole on the third sheet has disappeared due to a simi- 
lar mechanism as discussed in Ref. jSOll for the case of the 
A(1232); it is so far displaced from the physical axis that it has 
no visible effects any more. 

The conclusions from Ref. ifsoll concerning poles in T^^ ap- 
ply also for the Do partial wave: r'^^ can be large and non- 
perturbative, associated with dynamically generated poles in 
the complex plane. However, poles in T^^ are systematically 
displaced far in the complex plane once is added; only poles 
in the full T matrix have physical significance. As a result, the 
A^*(1520) appears as a clean resonance in Do, with no addi- 
tional structures in the amplitude. 



4. Summary and Conclusions 

The analytic properties of the nN Jiilich model have been de- 
termined. For this, a method for the analytic continuation for 
propagators with stable and unstable particles has been devel- 
oped. Through the deformation of the integration contour it is 
possible to investigate all unphysical sheets. A channel with 
stable particles induces one branch point and two sheets for the 
amplitude. For effective nnN channels, the pole of the unstable 
particle induces two additional branch points in the full propa- 
gator In total, every effective nnN channel induces three branch 
points and four sheets in the scattering amplitude. 

The pole positions and residues of the baryonic resonances 
up to a total spin of 7 = 3/2 have been extracted. Residues, bare 
and dressed couplings have been compared, showing that only 
the residue provides a well-defined expansion parameter free 
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of the model dependent decomposition into pole and non-pole 
part. For the D13 partial wave, poles in the non-pole term r'^^ 
are displaced far into the complex plane in the full amplitude 
r = -I- T^^; these findings are in line with the case of the f 33 
partial wave discussed in a previous study. 

Resonance interference was shown to play a crucial role in 
5 11. Taking this partial wave as an example, the quality of the 
extracted parameters has been shown; the residue terms alone 
provide already a good description of the full amplitude on the 
physical axis and also in the complex plane. 

The amplitudes of the Jiilich model are derived within a 
field theoretical approach from Lagrangians obeying chiral con- 
straints. Data are described to a high precision in the various 
partial waves. We claim that with these ingredients, in com- 
bination with a thorough treatment of the analytic properties, 
a reliable and precise extraction of pole positions and residues 
becomes possible. 
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A. The analytic structure of partial wave amplitudes 



Circular cut 



LH cut 




Figure 22: The analytic stmcture of the ttN partial wave amplitudes below 
threshold. Upper left: right-hand (RH). left-hand (LH), circular and short nu- 
cleon cuts in the j-plane. Upper right: Analytic structure at the intersection 
of circular and short nucleon (SN) cut. The short nucleon cut (SN) starting 
at follow on sheet 5 1 ends at bl^ip on sheet S2, whereas the short nucleon cut 
starting at fo2io„ on sheet S2 ends at hi^p on sheet S 1. Lower: The physical 
Riemann sheet is indicated with S 1 (blue surface), the analytic continuations 
along the circular and short nucleon cuts are indicated with S 2 (yellow surface) 
and S3 (orange surface), respectively. 



It is well-known from dispersion theoretical considerations 
1 62 , 3l that below the ttN threshold there are various additional 
cuts and branch points. These cuts lead to analytic continu- 
ations of the scattering amplitude; however, due to their rela- 
tively large distance to the physical scattering region, poles on 
those sheets are rather of academic interest and have little influ- 
ence on the partial waves. For completeness, in this Appendix 
the main features of the sub-threshold cuts are discussed. 

It is possible to identify the circular, short nucleon and other 
cuts with the partial wave projections of different diagrams of 
the present approach. First, there is the nucleon pole at z = 
niN in the f u partial wave, coming from nucleon s channel 
exchange. In the current approach, the bare nucleon mass and 
bare nNN coupling constant are renormalized by requiring the 
pole to be at mwand its residue to correspond to the physical 
ttA^A^ coupling III El- 

For u and f channel cuts, it is enough to consider the structure 
of M and f channel exchange processes contained in the present 
approach, projected to partial waves via Legendre polynoms P(. 
The two time orderings of the present TOFT formalism can be 
combined to structures appearing in the usual Feynman propa- 



gator according to 



Vu ^ f dx 

J M - mi, + le 

-I 

1 



t - m, + le 



(34) 



with u channel nucleon exchange and t channel exchange of 
particles with mass m, > 2m„ as given e.g. by the correlated 
In exchange in the cr and p channel, derived from A^A^ — > tttt 
pseudo-data via dispersion relations and using crossing sym- 
metry ||63, 64]. Consider the on-shell to on-shell nN nN 
kinematics as given by the transitions in Eq. ( l34l i. 

It is easy to see that y„ leads to the so-called short nucleon 
cut and both V,, and V, contribute to the so-called circular cut. 
The branch points of the short nucleon cut can be obtained by 
solving for the zero of the denominator of Vu at the borders of 
integration jc = + 1 , 



m|a-=±i = 0, 



(pi- 



(35) 
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where p„, p'^ are the final pion and initial nucleon momentum, 
respectively. The branch points of the short nucleon cut are 
situated above and below the z - m^. There are two short 
nucleon cuts as indicated in the upper right of Fig. |22]due to 
the presence of the circular cut that crosses the short nucleon 
cut. 

The projected potential V„ can be analytically continued 
along the circular cut (5 1 <-> 52), as indicated in the figure. 
The analytic continuation of the first Riemann sheet S 1 along 
the circular cut (CC), situated at 



keel = ^|'nl-ml (929 MeV) 
is indicated as 52. The two sheets are given by 



(36) 



Pf(x) 



{En - E„)- -2ql^(l - x) - ml 



y(52). r dx — ^ 

J (EN + E,)2-2q 



X) 

?on(l -x)-ml 



(37) 



with energies £, - ^qln + mj and the on-shell momentum qo„ 
from Eq. (|5]i. The short nucleon cut (SN) starting at blio„ on 
sheet S 1 ends at fo2up on sheet 52, whereas the short nucleon 
cut starting at b2\o„ on sheet S2 ends at bl^p on sheet 51. The 
branch points are situated at 



Z(^llow) = 


ml -ml 


(899 MeV, on SI) 


^ml + 2ml 


Z(^2up) = 


mN 


(939 MeV, on S2) 


z(b2iow) = 


ml - ml 


(919 MeV, on S2) 


mN 


Z(^lup) = 


J ml + 2ml 


(959 MeV, on SI). 



Note the branch point at z = mN on the sheet 52 is not on the 
physical sheet 5 1, where the nucleon pole is situated 0. 

The analytic continuation 5 1 — > 5 3 and 53^51 along the 
short nucleon cut is given by (only = is considered): 



,(Sl)/=0 m 



t=0 



^ for Im z > 



(39) 



for the part of the short nucleon cut with z < \zcc\- This con- 
tinuation is indicated with the orange surface in Fig. |22] The 
branch points of the short nucleon cut are logarithmic, i.e. in- 
finitely many sheets are connected at this branch point. See e.g. 
the surface 53 in Fig. |22]that connects to the next sheet (not 
drawn) below z < /jIiow We do not quote the corresponding 
analytic continuations, nor the continuation along the short nu- 
cleon cut 52 53. 



^^The rather involved structure of the short nucleon cuts discussed here is not 
obvious in Ref. d; in Ref. (el only one short nucleon cut is found. 



As Fig. |22] upper left, shows there is another short nucleon 
cut on the left-hand side obtained by z — > -z. There are addi- 
tional cuts induced by y„ from to ±; oo. Furthermore, there is 
an intermediate cut from -mN + m„ to mN - m„ coming from the 
kinematic factor k appearing in Eq. (l24l i. as well as a singularity 
at z = 0. Corresponding to the righthand, physical cut, there is 
the left-hand cut starting from z - -mN - m„ to -oo. These cuts 
do not come from the partial wave projection but arise from the 
unitarity contained in Eq. ([T]) and symmetry z — » -z. 

The t channel exchange V, from Eq. ( |34| ) contributes to the 
circular cut and the cut from to ±/oo. While « and t chan- 
nel exchange follow from crossing symmetry, there are many 
more interaction potentials contained in the present study. For 
example, the A(1232) u channel exchange leads to a short A cut 
situated between z - and the short nucleon cut. The poten- 
tials contained in the present study lfl9ll lead to cut structures 
and branch points beyond the structures discussed in this Ap- 
pendix. 

Additionally, the exchange potentials appear not only in on- 
shell kinematics, as discussed here, but also in half-off-shell 
and off-shell kinematics. The three body cuts induced by such 
kinematics have to be carefully taken into account in the pole 
search far from the physical axis as discussed at the end of Sec. 



In this Appendix we have shown that the general analytic 
structure of the pion-nucleon scattering amplitude Ii62t, i4i1 can be 
identified with different ingredients of the Jiilich model. In par- 
ticular, all cuts and branch points demanded by crossing sym- 
metry are present; circular and short nucleon cut, the interme- 
diate cut from -mN + m„ to mN - m„, the cuts from to +; oo, 
as well as the left-hand cut and a left-hand short nucleon cut 
obtained from z — » -z. 

These cuts and respective branch points can be indirectly vis- 
ible, e.g. in a sharp rise of the Re 5 1 1 amplitude close to the nN 
threshold as shown in Fig. [17] Indeed, most of the amplitude 
close to threshold is given by nucleon exchange plus correlated 
two-pion exchange, as we have tested. In models without these 
potentials, the sharp rise at the jiN threshold can lead to the 
appearance of sub-threshold poles that mimic the missing cut 
structure ||65|]. Such a behavior has been recently found in Ref. 
16611 , where the meson baryon interaction is solely given by the 
Weinberg-Tomozawa contact term. Such sub-threshold poles 
should not be confused with genuine states. 

References 

[1] nth Workshop on the Physics of Excited Nucleons, NSTAR 2007, Sep. 

5-8 2007, Eur. Phys. J. A 35, 253-331 (2008) 
[2] C. Amsler et al. [Particle Data Group], Phys. Lett. B 667, 1 (2008). 
[3] R. Koch and E. Pietarinen, Nucl. Phys. A 336, 331 (1980). 
[4] G. Hohler, Pion Nucleon Scattering, edited by H. Schopper, Landolt 

Bornstein, New Series, Group 9b, Vol. I (Springer, New York, 1983) 
[5] R. A. Arndt, W. J. Briscoe, I. I. Strakovsky and R. L. Workman, Phys. 

Rev C 74, 045205 (2006). 
[6] R. E. Cutkosky, C. R Forsyth, R. E. Hendrick and R. L. Kelly, Phys. Rev. 

D20, 2839 (1979). 
[7] R. L. Workman and R. A. Arndt, Phys. Rev. C 79, 038201 (2009). 
[8] A. V. Sarantsev, V. A. Nikonov, A. V. Anisovich, E. Klempt and 

U. Thoma, Eur. Phys. J. A 25, 441 (2005). 



21 



[9] A. V. Anisovich, A. Sarantsev, O. Bartholomy, E. Klempt, V. A. Nikonov 
and U. Thoma, Eur. Phys. J. A 25, 427 (2005). 
[10] D. M. Manley and E. M. Saleski. Phys. Rev. D 45, 4002 (1992). 
[11] T. Feuster and U. Mosel, Phys. Rev. C 58, 457 (1998). 
[12] G. Penner and U. Mosel, Phys. Rev. C 66, 05521 1 (2002). 
[13] V. Shklyar, H. Lenske, U, Mosel and G. Penner, Phys. Rev. C 71, 055206 

(2005) [En-atum-ibid. C 72, 019903 (2005)]. 

[14] T. P Vrana, S. A. Dytman and T. S. H. Lee, Phys. Rept. 328, 181 (2000). 
[15] S. Ceci, A. Svarc, B. Zauner, M. Manley and S. Capstick, Phys. Lett. B 

659, 228 (2008). 
[16] T. Sato and T. S. H. Lee, Phys. Rev. C 54, 2660 (1996). 
[17] Y. Surya and R Gross, Phys. Rev. C 53, 2422 (1996). 
[18] C. Schutz, J. Haidenbauer, J. Speth and J. W. Durso, Phys. Rev. C 57, 

1464(1998). 

[19] O. Krehl, C. Hanhart, S. Krewald and J. Speth, Phys. Rev. C 62, 025207 
(2000). 

[20] A. M. Gasparyan, J. Haidenbauer, C. Hanhart and J. Speth, Phys. Rev. C 
68, 045207 (2003). 

[21] B. Julia Di'az, T. S. Lee, A. Matsuyama and T. Sato, Phys. Rev. C 76, 
065201 (2007). 

[22] J. Durand, B. Julia Diaz, T. S. Lee, B. Saghai and T. Sato, Phys. Rev. C 
78, 025204 (2008). 

[23] H. Kamano, B. Julia Diaz, T. S. Lee, A. Matsuyama and T. Sato, Phys. 

Rev. C 79, 025206 (2009). 
[24] A. Dobado and J. R. Pelaez, Phys. Rev. D 47 (1993) 4883. 
[25] N. Kaiser, R B. Siegel and W. Weise, Phys. Lett. B 362 (1995) 23. 
[26] N. Kaiser, T. Waas and W. Weise, Nucl. Phys. A 612 (1997) 297. 
[27] J. A. Oiler and E. Oset, Nucl. Phys. A 620 (1997) 438 [Erratum-ibid. A 

652(1999) 407]. 
[28] E. Oset and A. Ramos, Nucl. Phys. A 635 (1998) 99. 
[29] J. A. Oiler, E. Oset and J. R. Pelaez, Phys. Rev. D 59 (1999) 074001 

[Erratum-ibid. D 60 (1999) 099906]. 
[30] U.-G. MeiBner and J. A. Oiler, Nucl. Phys. A 673, 311 (2000). 
[31] J. A. Oiler and U.-G. MeiBner, Phys. Lett. B 500 (2001) 263. 
[32] T. Inoue, E. Oset and M. J. Vicente Vacas, Phys. Rev. C 65 (2002) 035204. 
[33] D. Jido, J. A. Oiler, E. Oset, A. Ramos and U.-G. MeiBner, Nucl. Phys. A 

725, 181 (2003). 

[34] E. E. Kolomeitsev and M. R M. Lutz, Phys. Lett. B 585 (2004) 243. 
[35] S. Sarkar, E. Oset and M. J. Vicente Vacas, Nucl. Phys. A 750 (2005) 294. 
[36] M. Doring, E. Oset and D. Strottman, Phys. Rev. C 73, 045209 (2006). 
[37] S. Weinberg, Phys. Rev. 130, 776 (1963); 131, 440 (1963); 137 B672 
(1965). 

[38] V. Baru, J. Haidenbauer, C. Hanhart, Yu. Kalashnikova and A. E. Kudry- 
avtsev, Phys. Lett. B 586, 53 (2004); C. Hanhart, Eur. Phys. J. A 35, 271 
(2008). 

[39] N. Suzuki, T. Sato and T. S. Lee, Phys. Rev. C 79, 025205 (2009). 
[40] A. Mueller- Groeling, K. Holinde and J. Speth, Nucl. Phys. A 513, 557 
(1990). 

[41] J. Wess and B. Zumino, Phys. Rev. 163, 1727 (1967). 

[42] U.-G. MeiBner, Phys. Rept. 161, 213 (1988). 

[43] W.R. Prazier and A.W. Hendry, Phys. Rev. 134, B1307 (1964). 

[44] R. E. Cutkosky and S. Wang, Phys. Rev. D 42, 235 (1990). 

[45] D. Gamermann and E. Oset, Eur. Phys. J. A 33, 1 19 (2007). 

[46] M. Albaladejo and J. A. Oiler, Phys. Rev. Lett. 101, 252002 (2008). 

[47] T. Sato and T. S. Lee, arXiv:0902.3653 [nucl-th]. 

[48] S. Ceci, A. Svarc and B. Zauner, Eur. Phys. J. C 58 (2008) 47. 

[49] S. Capstick et al, Eur. Phys. J. A 35, 253 (2008). 

[50] M. Doring, C. Hanhart, R Huang, S. Krewald and U.-G. MeiBner, 

prXiv:0903.1781 [nucl-th]. 
[51] I. R. Afnan and B. Blankleider, Phys. Rev. C 22, 1638 (1980). 
[52] A. Matsuyama, T. Sato and T. S. Lee, Phys. Rept. 439, 193 (2007). 
[53] M. Haftel, R Tabakin, Nucl. Phys. A 158, 1 (1970). 
[54] 1. Caprini, G. Colangelo and H. Leutwyler, Phys. Rev. Lett. 96, 132001 

(2006) . 

[55] V. K. Magas, E. Oset and A. Ramos, Phys. Rev. Lett. 95, 052301 (2005). 
[56] G. Hohler, nN Newsletter 9 1 (1993). 

[57] R. Arndt, W. Briscoe, 1. Strakovsky and R. Workman, Eur. Phys. J. A 35, 
311 (2008). 

[58] M. Dugger et al. [CLAS Collaboration], Phys. Rev. C 79, 065206 (2009). 
[59] R. A. Arndt, W. J. Briscoe, 1. 1. Strakovsky, R. L. Workman and M. M. Pa- 
van, Phys. Rev. C 69, 035213 (2004). 



[60] M. Doring et al, in preparation. 

[61] S. Sarkar, E. Oset and M. J. Vicente Vacas, Eur. Phys. J. A 24, 287 (2005). 
[62] S. C. Frautschi, J. D. Walecka, Phys. Rev. 120, 1486 (1960). 
[63] C. Schiitz, J. W. Durso, K. Holinde and J. Speth, Phys. Rev. C 49, 2671 
(1994). 

[64] C. Schiitz, K. Holinde, J. Speth, B. C. Pearce and J. W. Durso, Phys. Rev. 

C51, 1374(1995). 
[65] R. A. Arndt, private communication. 
[66] M. Doring and K. Nakavama, larXiv:0906.2949] [nucl-th]. 



22 



